QCD radiative corrections to Xcj production at hadron colliders 
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To clarify the outstanding problem in charmonium production that existing theories cannot ex- 
plain the observed cross sections of Xcj{J = 0, 1, 2) and ratio R Xc = o- Xc2 /a Xcl w 0.75 (in contrast to 
the spin counting value 5/3) at the Tevatron, we study the complete next-to-leading order radiative 
corrections in nonrelativistic QCD, and find next-to-leading order contributions of 

3p[l] 

are more 

important than leading order at high pr, and decreases slower than , implying a natural 
explanation for the R Xc puzzle. By fitting R Xc , the predicted cross sections of Xcj are found to agree 
with data. The result indicates color-octet contribution is crucially needed, thus providing a unique 
test for heavy quarkonium production mechanisms. Feed-down contributions of Xcj to prompt J/ip 
production are estimated to be substantial, about 30 — 40% at pr — 20 GeV. Production of Xcj 
(J=0,l,2) at the LHC is also predicted. 
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Heavy quarkonium production remains a challenging 
problem in understanding quantum chromodynamics. 
Among others, the puzzle of J/ip(ip') production cross 
sections and polarizations at large transverse momentum 
Pt at the Tevatron is crucial in testing the color-octet 
(CO) and color-singlct (CS) mechanisms in NRQCD 
and other mechanisms (for reviews see e.g.Q). The P- 
wave charmonia Xcj production is equally important, 
since they give substantial feed-down contributions to 
the prompt J/ip production through decays Xcj — > iJ/^P- 
An even more important issue in Xcj production concerns 
the ratio of production rates of Xc2 to xa- 

The CDF Collaboration measured the ratio 
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to be about 0.75 for pt > 4 GeV, and even smaller 
whenpT becomes larger^). But at leading order (LO) in 
a s , NRQCD predicts the Xc production cross sections to 
>[i] 



LO predictions Q for J/ip exclusive and inclusive produc- 
tion in e + e~ annihilation at B factories can be resolved 
at NLO within CS channels 0. For J/ip photoproduc- 
tion at HERA, complete calculations^) slightly favor the 
presence of CO contributions. For J/ip production at the 
Tevatron, the NLO correction in CS channels enhances 
the cross section at large px by 2 orders of magnitude 
[|9, 10 and reduces the discrepancies between theory and 
experiment. So, a crucial issue in charmonium produc- 
tion is whether the CO contributions are still needed. To 
clarify this, we must go beyond S-wave and study P-wave 
quarkonia, e.g., Xcj production. 

In view of the urgency, we study the NLO QCD cor- 
rections to XcJ hadroproduction in this work. Within 
NRQCD factorization, the inclusive cross section for the 
direct Xcj production in hadron-hadron collisions reads 



da[pp Xcj + X] = da[(cc) n ] 
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scale as l/p T in the CS "T) ' channels yet scale as l/p T in = J2 I & x i& x 2G l / p G ] / p x da[i + j ^ (cc) n + X) (O. 
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the CO S[ 1 channel. Thus S[ would dominate at large 

Pt-, predicting the ratio to be 5/3 by spin counting^ 
which is much larger than the measured value 0.75. 
Meanwhile, the color-evaporation model (CEM) predicts 
the ratio to be always 5/3 in all orders of a s simply 
based on spin counting. It seems as if no existing theory 
agrees with the measured R Xc - However, there could be 
a good chance for NRQCD to explain this puzzle, be- 
cause the next-to- leading order (NLO) contributions in 
a s will change the large px behavior of cross sections. 
In particular, contributions of CS T'jp channels scale as 
l/p T at NLO, more important than l/p T at LO. So it 
is necessary to study Xcj production at NLO to see how 
the value of R Xc can be changed at large px- 

Another issue concerns the CO mechanism, which has 
been studied extensively in S-wave charmonium J/ip pro- 
duction. Large discrepancies between experiments and 
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where p is either a proton or an antiproton, the indices 
i,j run over all the partonic species and n denotes the 
color, spin and angular momentum (L n ) of the interme- 
diate cc pair. In this work, we calculate the cross sec- 
tions at NLO in a s and LO in v (the relative velocity 
of quark and antiquark) . So only two channels and 
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S[ in the present calculation are involved. The long- 
distance-matrix elements (LDMEs) (C^ c J ) are related to 
the transition probabilities from the intermediate state 
(cc)„ into Xcj, and are governed by the nonperturbativc 
QCD dynamics. Note that our definition of CS LDME 
{O x ° J (?P [ j ] )) is different from that in Rcf.ft| by a factor 

of l/(2iV c ). And (O x " J ( 3 Pj 1] )) (J=0,l,2) are related to 
just one matrix element by spin symmetry. 

We use FeynArts]12|| to generate Feynman diagrams 
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FIG. 1: Representative Feynman diagrams for \cj hadropro- 
duction at LO and NLO. The gluon-gluon and gluon-quark 
subprocesses are all included. In (b)-(f) both color-singlet 
3 -Pl and color-octet 3 £>[ 8 ' channels contribute. 



and amplitudes. Some representative diagrams are 
shown in Fig. |l|. There are generally ultraviolet (UV), 
infrared (IR) and Coulomb singularities. Conventional 
dimensional regularization (CDR) with D = 4 — 2e is 
adopted to regularize them. We have checked analyti- 
cally that all singularities are canceled exactly. 

The UV divergences from self-energy and triangle di- 
agrams are removed by renormalization. The renormal- 
ization constants Z m , Z2, Z21 and Z3, which correspond, 
respectively, to the charm quark mass to c , charm- field 
ipe, light quark field ip q and gluon field A® are defined in 
the on-mass-shcll(OS) scheme, while Z g corresponding to 
the QCD gauge coupling a s is defined in the modificd- 
minimal-subtraction(MS) scheme: 
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where N f 



e) is a overall factor in our 
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calculation, /3rj(n/) = ^-Ca ~ |Tj?n/ is the one-loop co- 
efficient of the QCD beta function, nf = 3 is the number 
of active quark flavors and fi r is the renormalization scale. 

IR singularities arising from loop integration and phase 
space integration of the real correction partially can- 
cel each other. The remaining singularities of the S- 
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wave state S[ are absorbed by the proton parton- 

distribution- functions (PDFs), while that of 3 Pj 1 ' states 

are absorbed by both the proton PDFs and (O x "°( 3 sf ] )). 
We extract poles in the real corrections using the phase 



space slicing method jl3|. We note that to correctly get 
the soft poles, the eikonal current should be taken before 
the heavy quarks are coupled to states with definite quan- 
tum numbers. We verify that our results are independent 
of the two cuts introduced in the phase space slicing for 
each subgroup associated with a real correction process. 

There are thousands of diagrams which arc handled by 
our self-written Mathematica program and then changed 
into C++ code to perform convolution and phase space 
integration. To perform the calculation, two different 
methods are used, resulting in two independent computer 
codes. In one of our methods, the virtual corrections are 
calculated numerically and using QCDLoopjli] to separate 
the singularities and finite result, while the real correc- 
tions are calculated using the helicity method. In the 
second method, we expand the singularities analytically 
and use LoopTools|l5j to get the finite result, while the 
real corrections are calculated by directly squaring the 
amplitudes. Results of the two methods are found to 
coincide with each other. 

Essentially, in our calculation the most complicated 
part is the loop-correction for P-wavc channels, where the 
derivation of the amplitude with respect to the relative 
momentum of heavy quarks Q and Q is involved. This 
derivation is equivalent to having one additional massive 
vector boson in the calculation, which causes complicated 
tensor loop integrals and entangled pattern of infrared 
(IR) singularities. Based on the formula in Ref. [IlJ], wc 
developed a new method to perform NLO corrections for 
processes involving bound states. This method makes the 
calculation of P-wave heavy quarkonium hardroproduc- 
tion at NLO become possible. The details of our method 
will be presented elsewhere Jl6|. 

Thanks to the LHAPDF interface |L7), the CTEQ6L1 
and CTEQ6M PDFs are used for LO and NLO calcula- 
tions respectively. The charm quark mass is set to be 
m c = 1.5 GeV, while the renormalization, factorization, 
and NRQCD scales are /i r = /// = tot and /i\ = m c , 
where tot = ^ + 4m^ is the XcJ transverse mass. The 
center-of-mass energies arc 1.96 TeV at the Tevatron 
RUN II and 14 TeV at the LHC. To estimate theoreti- 
cal uncertainties, we vary /i r and [if from tjit/2 to 2tot 
and choose m c = 1.5 ± 0.1 GeV. Note that there exists 
a large cancellation for errors between fitting data and 
predictions. Thus to avoid double counting, our strategy 
is to consider all errors in the fit procedure and include 
them in fitted parameters, while when predicting new 
quantities, we just consider errors due to fitted parame- 
ters. 

We begin our numerical analysis by comparing the 
short-distance coefficients d(i[ 3 P^] of the ^j 1 ' channels 

and dafsf ] ] of the 3 sf channel at NLO. Wc find that at 
NLO the most significant change is the large px behav- 
ior in the Tj channels: it scales as 1 /p T at large pr at 
NLO, in contrast to l/p T at LO. As shown in Fig. ||, at 
LO the Pj channels arc negligible at high px, whereas 
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FIG. 2: The ratio of dafp]^ ] ] to A&fsf ] ] at NLO as a function 
of pr at the Tevatron. The cut \y Xc , \ < 1 is chosen to compare 
with the CDF data of Ref. §. 



1.0 




0.0 1 

5 10 15 20 25 30 

Pr of J/ip (GeV) 



FIG. 3: Transverse momentum distribution of ratio 
R Xt .j 'Rj/^~/ at the Tevatron with cut |y Xc jl < 1- The CDF 
data is taken from Ref. |J. The lower and upper bounds of 
LO and NLO are constrained by 0.24 < r < 0.33. 



at NLO the channels are comparable to the 
channel even at pr ~ 50 GeV. Another important but 
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unique feature is that the Pj channels give large values 
at highpT, but surprisingly with negative signs. We note 
that the negative values are not caused by the choice of 
//a > though this may affect their absolute values. In fact, 
detailed investigation reveals that the negative values are 
originated from the renormalization scheme (RS) for the 
NRQCD LDMEs (C Xrf) ( 3 s[ 8] )) . The RS in this work is 
the conventional MS scheme. One may use another RS 
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to get different values of da [ Pj ] , but this should not 
change the physical result, because the RS dependence 
of short-distance coefficients of ^Pj 1 ' are canceled by the 

RS dependence of (O Xc0 ( S[ J )), and the final physical 
results are independent of the choice of it. Especially, 
dc[Xcj] an d R Xc ar c independent of the RS and ^a, and 
their values should always be positive. 

From Fig. |^ we sec that the ^p] 1 ' channel decreases 
slower than the P^ channel. Considering also that Pj 
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channels are comparable to the S[ channel at high px, 
we may naturally explain the CDF data[|| that the pro- 
duction rate of xa exceeds that of \c2 at high pt- We 
define the ratio 

r - <Q**fo W )> I (4) 

The bound of r > 0.24 is needed to get a positive pro- 
duction rate of Xco at high pr, as shown in Fig. |2[ With 
this definition, we can give an asymptotic expression of 
R Xr at the Tevatron with |y Xc / | < 1 : 

5 rdal 3 S [ * ] ]+daM 1] ] 5 r - 0.20 
Xc 3 r ^] ]+ ^W] 3r-0.16' W 

where the numbers —0.20 and —0.16 can be read from 
Fig. ||. Because of these two numbers, R Xc must be 
smaller than 5/3 at high pr- By fitting the data [||, 
we get r = 0.27^°;°^° °^°;^ w 0.27 ± 0.06, where the 
errors come from data, scale dependence and m c depen- 
dence respectively. In the fitting procedure, the mass 
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FIG. 4: The same as Fig. | but for LHC with cut \y Xc J | < 3. 

difference between Xcj and J/tp is also considered. It 
can be found that the value of r is compatible with the 
naive velocity scaling rule r ps 0(1) 0. 

As shown in Fig. ||, the NLO predictions present a 
much better compatibility with data than LO, where r 
is constrained by 0.24 < r < 0.33 and Rj/^y = B(xd —> 
J/ipj)/B(xc2 -> J/il>l) = 1-91 as in Ref.@. In Fig. |, 
we show predictions for R Xa at the LHC. Note that, the 
ratio R x is sensitive to r at NLO. Thus measurement of 

Ac 

R Xc at high pr will give a strict constraint on r. 
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FIG. 5: Transverse momentum distribution of Xcj feed down 
to J/ip at the Tevatron RUN I and LHC, where RUN I data 
is taken from Ref. |l9| . 
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To predict the production cross sections, we still need 
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a CS matrix element (O Xc0 (Pj )). As a widely adopted 
choice, we take the potential model result |i?p(0)| 2 = 
0.075 GeV 5 (see the B-T model in Q) as our input pa- 
rameter. Then we compare our prediction for \ c j pro- 
duction at the Tevatron with RUN I data in Fig. |5|, where 
= 1.8 TcV and \y XcJ \ < 0.6. Because of a large 
negative correction for CS channel at large pt, NLO re- 
sults decrease faster than LO, and give a more reason- 
able interpretation for experimental data. The differ- 
ences between \cj production at the Tevatron RUN II 
with VS = 1.96 TeV and that of RUN I are within a 
factor of 20 percent for all pr region, so we do not present 
them here. At the Tevatron, the feed-down contribution 
from Xcj has a rather large proportion of the total J/ ip 
production cross section. If we choose r = 0.24(0.27), 
the proportion is 25(30)% at pr = 10 GeV, and reaches 
to about 30(40)% at pr = 20 GeV. Because of this large 
proportion, \cJ feed-down will have an important ef- 
fect on J /if; polarization in J '/tp prompt production, and 
should be further clarified when dealing with the J/ip po- 
larization puzzle. The prediction for XcJ production at 
the LHC is also presented in Fig. [| where VS = 14 TcV 

an d \y X cj I < 3. 

In summary, we calculate NLO QCD corrections to 
XcJ production at the Tevatron and LHC, including both 
CS and CO channels. We find ^pj 1 ' channels give large 

3 [11 

contributions at high pr and P{ decreases slower than 
^Rp' ; then the measured ratio of R Xc at the Tevatron can 



be naturally explained. Moreover, our result shows that 
the measured R Xc disfavors CEM, but favors NRQCD. 
By fitting the observed ratio R Xc , we extract the ratio of 
CO to CS matrix element to be r = 0.271q'q3, which is 
used to predict the production rates of Xcj at the Teva- 
tron RUN I and leads to a good agreement with data. As 
a result, for the first time, the observed rates of XcJ and 
ratio R Xo are explained simultaneously. Our result may 
also be used to predict the prompt J/ip production and 
shed light on the J /ip polarization puzzle. 

We emphasize that in Xcj production the NLO cor- 
rections already scale as l/pj., which is the leading pt 
behavior, and the NNLO and other corrections are no 
longer important, because they are suppressed either by 
a s or by v 2 relative to NLO contributions. This differs 
markedly from J /ip production, where NNLO corrections 
could be even more important than NLO at large pr [^0) ■ 
As a result, we have picked up the most important con- 
tributions in Xcj production, and given a good descrip- 
tion for this process. This work gives strong support to 
the color-octet mechanism, and provides further tests for 
NRQCD and heavy quakonium production mechanisms 
at the LHC. 
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